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We study the linewidth of the surface plasmon resonance in the optical absorption spectrum of 
metallic nanoparticles, when the decay into electron-hole pairs is the dominant channel. Within 
a semiclassical approach, we find that the electron-hole density-density correlation oscillates as a 
function of the size of the particles, leading to oscillations of the linewidth. This result is confirmed 
numerically for alkali and noble metal particles. While the linewidth can increase strongly, the 
oscillations persist when the particles are embedded in a matrix. 
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I. INTRODUCTION 

The lifetime of collective excitations in finite quan- 
tum systems is a fundamental problem spamiing from the 
physics of nuclei to that ofrfnetallic clustersEJ. Giant res- 
onances, in the former caseB, and surface plasmon, in the 
latterEI, decay by different mechanisms according to the 
physical regime of the excitation. For large metallic clus- 
ters (with radius i? ^ 10 nm) radiation daruping is the 
main limiting factor to the plasmon lifetimcQ. At small 
sizes (0.5 nm ^ i? < 2nm) the decay into electron-hole 
pairs (Landau damping) dominates, while for intermedi- 
ate sizes both effects compete. For very small clusters 
(i? ^ 0.5jim) a strong temperature dependence of the 
linewidthd suggests that the coupling of the collective 
excitation to thermal vibrations of the atomic lattice is 
an important effect. 

Except for the case of very small clusters, where the 
atomic arrangements play a crucial role, the theoretical 
analysis of the surface plasmon is usually based on the 
jellium model, where conduction electrons move in a uni- 
formly charged background. Within such a framework, 
the Landau damping occurs when there are particle-hole 
excitations with an energy close to that of the collective 
excitation. Therefore, this mechanism is strongly depen- 
dent on the density of states of the conduction electrons 
confined in the cluster. 

Usiflg a linear response approach, Kawabata and 
KuboQ proposed a total linewidth F = Fi-|- AF(i?), where 
the intrinsic width Fi is given by bulk properties and 
the size-dependent contribution AF(i?) is proportional 
to 1/ R. Such an R dependence would also result from 
the identification of the diameter of the particle with 
the electron mean-free-path. However, as stressed in the 
original work, this identification is not completely cor- 
rect, since the quantized states exist due to the confine- 
ment Mid cannot be considered to be scattered at the 
surfacetllj. On the other hand, for interacting electrons 
confined in a parabolic potential, Kohn's thjeoremCI im- 
plies an infinitely-lived collective excitationB, such that 



the linewidth of the plasmon can be thought as a measure 
of the failure in separating the collective (center-of-mass) 
and relative coordinates for the electron system. 

Further improvements in the-, linear response the- 
ory by Yannouleas and BrogliaB and by Barma and 
SubrahmanyamllJ recovered the \/R dependence of 
AF(i?) in thej-asymptotic limit of large ii— . Numerical 
calculationsEM and existing experimentsaO on free in- 
termediate-size clusters are in good agreement with the 
values of AF(i?) yielded by those calculations. For small 
clusters the emergence of a non-monotonous size-depend- 
ent structure in the lineshape (.in. exuCTmyntalj'Efllij as 
well as in numerical calculationst2ltiHl3ElHl3) has_heen at- 
tributed to shell effects (see also Fig. 10 of Ref.Eja). We 
provide in this work a detailed theoretical understand- 
ing of these effects and determine analytically the size- 
dependence of the plasmon lifetime for free and embed- 
ded clusters. 



Using a semiclassical approach, we explain how a 
subtle consequence of electronic shell effects, namely 
angular-momentum dependent electron-hole density-den- 
sity correlations, leads to an oscillatory size- dependence 
of the plasmon linewidth due to Landau damping in small 
alkaline metal clusters. This result is confirmed by nu- 
merical calculations using the iiipe-dependent local den- 
sity approximation (TDLDA)EilJij. We show that our 
results are relevant also for noble metals and clusters em- 
bedded in matrices. Nanoparticles embedded in a matrix 
are interesting from the experimental point of view since 
the position of the resonance can be tuned by changing 
the dielectric surrounding. They constitute a compli- 
cated system, where the plasmon lifetime strongly de- 
pends on the properties of the matrix, aiijiffect that has 
been called chemical interface dampingljEd. The appli- 
cation of the TDLDA, within a simple model of the po- 
tential at the surface, confirms that the size-oscillation of 
the level-width predicted by our analytical method per- 
sists even for embedded clusters. 



II. LANDAU DAMPING IN NANOPARTICLES 

The nanoparticles we consider are much smaher than 
the wavelength of the radiation, and the optical absorp- 
tion cross-section at frequency u is given by 
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where the dipole matrix element is taken between the 
electronic many-body ground state |0) and the excited 
states \j) with energy Ej (the z-axis is chosen parallel 
to the electric field direction). The spectrum is domi- 
nated by the collective surface plasmon excitation which 
is characterized by an oscillation of the center of mass of 
the electrons with respect to the.lattice of positive ions. 
Following Mie's classical theoryQ, for a nanoparticle in 
vacuum one expects the surface plasmon resonance at 
the frequency ojm = Wp/-\/3, where Wp is the bulk plas- 
mon frequency of the material. 

When the confinement is not parabolic the surface 
plasmon is not an eigenstate of the many-body Hamil- 
tonian and its coupling to the particle-hole excitations 
leads to a broadening of the corresponding resonance in 
the spectrum (see the inset of Fig. P). 
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A common approach to the plasmon lifetime consists 
in treating the collective excitation as an external pertur- 
bation which can give rise to the creation of electron-hole 
excitationsEl. Within this framework, Fermi's Golden 
Rule yields the linewidth 



AT{R) = 2iTY,\(p\SV\h)f 6{hujM 
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where \p) and \h) are electron and hole states in the 
self-consistent field, respectively, and 5V is the dipole 
field due to the surface plasmon. More rigorous ap- 
proaches, like discrete-matrix RPA would lead after cer- 
tain approximationsti, to an equivalent separation of col- 
lective and particle-hole states. 

In the case of spherical symmetry we can work in the ef- 
fective one-dimensional problem of the radial coordinate 
for each value (in units of h) of the angular momentum 
L < p ^ Ry/2ME/h, E and M being the energy and 
the effective mass of an electron, respectively, and R de- 
noting the radius of the nanoparticle. Assuming that the 
confinement and the interactions lead to hard walls at 
radius R in the self-consistent field, |and integrating over 
the electron-hole states, one obtainsEl 



Ar(i?) = C7 f dEj2 Y. 

i L L'=L±l 



(2L + 1){2L' + 1) < L, 0; 1, 0|L', >2 E{E - hujM)dL{E)dL'iE ~ Hum) , (3) 



where < L, 0; 1, 0\L', > is a Clebsch-Gordan coefficient, 
7 = {27rh^)/{3NM^u;MR'^), C = AMR'^/h^, and dLiE) 
is the one-dimensional density of states. 



III. SEMICLASSICAL THEORY 

Our approach consists in the evaluation of (g|) by using 
the semiclassical expressiontil 
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which decomposes dL{E) in its smooth and oscillating 
components. Here, t{E,L) = 'h\/ p^ — LP- jE is the pe- 
riod of the one-dimensional motion, 



S{E, L) = 2n yjp^ - i2 _ L arccos(L/p) 



is the corresponding action, and the sum runs over all 
numbers of repetitions r. Integration over i, Poisson's 



summation rule, and a stationary phase apppapcimation, 
lead to the semiclassical Berry- Tabor formulac3, express- 
ing the density of states of a sphere in terms of classi- 
cal periodic trajectories. Taking only into account the 
smooth part of dL{E) in Eq. (g), we recover the well- 
known 1/ R dependences 
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where ^ = htOM/E-p , g{S^) is a smoothly decreasing func- 
tion with g(0) — 1, and vp (E-p) is the Fermi velocity 
(energy). 

Most interestingly, the oscillating part of the density 
of states gives rise to additional new features, namely an 
oscillation of AT{R) around To{R) as a function of the 
radius. While the cross products between the smooth and 
the oscillating parts of d^^E) and dLi{E — Tiujm) can be 
neglected after the energy integration, the product of the 
two oscillatory components, after doing the L summation 
with the aid of Poisson's summation rule, yields 
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with p' = ^ p^ — (fcpi?)^^. The hmits of the p-integral are given by p^m — ixia.x{kpR, kpR^) and pmax = kpR^/l + ^. 
We have defined /(y') = y' + 1 for y' = y — f and /(y') = y' in the case y' — y +1 (fcp is the Fermi wave- vector) . The 
dominant terms correspond to m = 0, r = — r' and y ^ y' + 1 = p/{p — p'), yielding 
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Since p is of the order of kpR, and kpR^ ^ 1, the argu- 
ment of the cosine is close to 2rkpR(,, and (|g) yields a 
contribution to AF which oscillates as a function of the 
radius of the nanoparticle. It is important to note that 
this oscillatory behavior is due to the angular-momentum 
restricted electron-hole density-density correlation, and 
therefore more subtle than the well-known oscillations 
of the density of states with the energy. 

Fig. m shows the result for AF = Fq H- Fqsc which 
is obtained by numerically evaluating the p-integral in 
Eq. (^ . The typical amplitude of the oscillations can be 
estimated by the value F^'^^'^ ^ Gy^h^ / {M Rl^ y^k^R^) , 
valid in the limit kpR^,^ ^ 1- This shows that the oscil- 
lating term represents an important correction to Tq{R) 
for small particle sizes. 

Higher repetitions (r > 1) of the periodic orbits are 
suppressed in Eq. ^) because of the prefactor l/^/r, and 
the oscillations of the cosine as a function of p. Includ- 
ing a small intrinsic dccoherence rate Fj accounts for the 
coupling to the cnviroiwHcnt and leads to an additional 
factor exp{—T[rR/ hvFf3, which further reduces the os- 
cillatory contributions from long trajectories with large 
repetition numbers r. 




FIG. 1: Linewidth Ar(7?), as a function of the radius (in 
units of the Bohr radius ab ~ 0.53A), for Na nanoparticles 
containing between 20 and 832 atoms. The dots correspond- 
ing to TDLDA calculations are compared to the result of 
Eq. (|) (full line). The smooth term To{R) is given by Eq. (|) 
(dashed). Inset: Logarithm of the numerically calculated ab- 
sorption cross-section for Naisg, showing a pronounced sur- 
face plasmon resonance, fitted by a Lorentzian (dotted line). 
The excited states are indicated by tick marks and their oscil- 
lation strengths are given by the height of the vertical lines. 



IV. NUMERICAL CALCULATIONS 

In order to study the relevance of the oscillatory be- 
havior in more general situations, we performed quan- 
tum numerical calculations within a spherically syrnrnet.- 
ric jellium model, using the linear response TDLDAEjo. 
The diagonalization of the Hamiltonian in the space of 
particle-hole excitations amounts to use the Random 
Phase Approximation in the self-consistent field with a 
local exchange potential. Further refinements of the bh- 
merical approach (like the Self-Interaction CorrectiorO) 
do not provide substantial changes in the effects we are 
interested in. Because of the spherical symmetry of the 
model, we considered nanoparticle sizes corresponding to 
magic numbers of atoms where our numerical approach 



Ji. 



is expected to work the bes 



A. Alkaline metals 

A typical spectrum of the photo absorption cross- 
section (Eq. (|l|)) is shown in the inset of Fig. |^ for the ex- 
ample of a sodium cluster, where a non-zero Fi smears the 
singularities of the spectrum. The collective excitation 
can be identified as the main peak of the spectrum (note 
that the logarithm of a is plotted). The vertical lines 
are essentially given by single electron-hole excitations, 
except for the one at the center of the surface plasmon 
resonance, which is the superposition of many such ex- 




FIG. 2: The dots represent the numericaUy obtained surface 
plasmon hnewidth for the case of silver nanoparticles embed- 
ded in an argon matrix, as a function of their radius R in units 
of the Bohr radius. The dashed and the solid lines represent 
ro(-R) and Ar(i?), respectively, reduced by a global factor of 
3 with respect to Eqs. (W) and (pi) (see text). 



ical results. This discrepancy is not surprising given the 
stringent approximations involved in both results, which 
are expected to be less applicable for embedded noble 
metal clusters, than for free alkaline metals. Neverthe- 
less, the analytically calculated rosc(-R) exhibits oscilla- 
tions with the correct period, and its relative size (with 
respect to Tq{R)) agrees with the numerical results. 

Reactive (non-inert) matrices like glass are often stud- 
ied experimcntallyll3, but still more difficult to analyze. 
By taking intekiaccount a high-energy conduction band in 
the insulatpiri, wc extended our numerical calculations 
to this cascEJ. These electronic states provide additional 
decay channels for the surface plasmon and its lifetime is 
considerably reduced. Our results for the Landau damp- 
ing contribution to the linewidth of Ag nanoparticles em- 
bedded in Si02 amount to about one half of the exper- 
imental linewidth, but still show strong size-dependent 
oscillations, as in the previous cases. 



citations. The surface plasmon linewidth F, represented 
by thick dots in Fig. ^, is obtained from a Lorentzian fit 
to the main peak. We chose an intrinsic width Fj larger 
than the typical distance in the electron-hole spectrum 
and verified that AF(i?) is approximately independent 
of the precise value of Fj. Confirming our analytical cal- 
culations, at small R the numerical results show strong 
deviations from Fo(i?) (dotted line) and a non-monotonic 
behavior as a function of the radius. Similar results are 
obtained for other alkaline metals. It should be men- 
tioned that the phase of the semiclassical result from 
Eq. (|6|) is meaningful up to a constant shift which was 
adjusted according to the numerical data. 



B. Noble metals and embedded clusters 

Noble metals are more difficult to describe, due to 
band-structure effects and the fact that the comparison 
with experiments has to take into account the confine- 
ment of the particles in a matrix. In our numerical ap- 
proach we included the contribution ed of the d-electrons 
to the dielectric function of the metal particles, as well as 
the dielectric constant Cm of the matrix. For the result- 
ing residual interaction we used a multipolar expansion 
as in Ref.EZl. We have checked that iacluding an intersti- 
tial spacing (the so-called 3-e modelc3) does not change 
significantly our results. 

In Fig. g we show the numerical results (thick dots) 
for silver nanoparticles embedded in an inert matrix 
(Ar, Em « 1.75), exhibiting pronounced oscillations of 
AF(i?) with the particle size, much like in the case 
of alkaline metals. The smooth part of the linewidth 
Fo(i?) (Eq. (|j)), and our semiclassical result for Fosc(-R) 
(Eq. (0)) are an overall factor of 3 larger than the numer- 



CONCLUSIONS 



In conclusion, we have found a mechanism leading 
to an important oscillatory contribution to the Landau 
damping linewidth of the collective surface plasmon ex- 
citation in small metallic nanoparticles. The oscillations 
of the linewidth as a function of the particle size arise 
from the oscillations of the angular-momentum restricted 
electron- hole density-density correlations, that we cal- 
culate within a semiclassical approach. We obtain a 
quantitative agreement between analytical and numeri- 
cal (TDLDA) results for free alkaline metallic clusters, 
and a qualitative agreement for more complicated cases. 

A systematic experimental investigation of the line- 
width oscillations exten d ima . th e results for specific situ- 
ations reported in Refs.E3iijii3 would provide an impor- 
tant step towards the understanding of the rich physics 
involved in the collective excitations of finite systems. A 
small inhomogeneous broadening will not kill the effect 
if the size dispersion of the nanoparticles is smaller than 
the period of the oscillations (approximately {k-p£)~'^ ^ ac- 
cording to Eq. (ph). Moreover, recent experimental de- 
velopments make it possible to study the homogeneous 
broadening of th p s ur face plasmon also in electrically un- 
charged clusterstilEHl. In the case of clusters embedded 
in a matrix, the oscillations of the width of the plasmon 
can be seen not only as a function of the size but also 
as a function of the dielectric constant of the environ- 
ment. Changing em shifts the position of the plasmon, 
without considerably modifying the single-particle part 
of the spectrum. Therefore the plasmon linewidth varies 
according to the position of the resonance with respect 
to the particle-hole excitation spectrum. 
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